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Colloidal crystal structures with complexity and diversity rivaling atomic and molecular crystals 
have been predicted and obtained for hard particles by entropy maximization. However, so far 
homochiral colloidal crystals, which are candidates for photonic metamaterials, are absent. Using 
Monte Carlo simulations we show that chiral polyhedra exhibiting weak directional entropic forces 
self-assemble either an achiral crystal or a chiral crystal with limited control over the crystal handed¬ 
ness. Building blocks with stronger faceting exhibit higher selectivity and assemble a chiral crystal 
with handedness uniquely determined by the particle chirality. Tuning the strength of directional 
entropic forces by means of particle rounding or the use of depletants allows for reconfiguration 
between achiral and homochiral crystals. We rationalize our findings by quantifying the chiral¬ 
ity strength of each particle, both from particle geometry and potential of mean force and torque 
diagrams. 


Controlling -and understanding the origins of- chiral¬ 
ity is a major goal in the physical and chemical sciences. 
It has longstanding implications that range from the de¬ 
sign of effective pharmaceuticals to explaining biological 
homochirality and the beginnings of life In recent 
years the expectation that chiral materials can provide 
a route for photonic metamaterials [5] has renewed the 
interest in chirality, now from a materials design perspec¬ 
tive. A key remaining challenge for the development of 
metamaterials [3] is the ability to use scalable techniques 
to create crystals with unique chirality. Of particular 
importance are enantioselective processes, which select a 
pre-defined handedness of the crystal, thereby maximiz¬ 
ing optical activity. 

While several routes have been reported to assemble 
structures with chiral order, from exploitation of elec¬ 
trostatic H], anisotropic [5] and isotropic [5] interac¬ 
tions, packing CHS], spontaneous symmetry breaking in 
twisted ribbons [ini E], to self-sorting in liquids [13] . 
the lack of chiral forces in these examples preclude ho¬ 
mochirality. Even if chiral forces such as circularly polar¬ 
ized light m are present, or when assembling naturally 
handed viruses bacteria m or nanoparticles with 
pre-sorted chiral coatings | 16 j . or even by judiciously 
designing chiral building blocks [T7H20|, the final chiral 
structures have, so far, been restricted to liquid crys¬ 
talline and low-dimensional arrangements. Growth of 
three-dimensional homochiral colloidal crystals has yet 
to be demonstrated. In fact, the only examples of chiral 
crystals with realized applications as polarization sensi¬ 
tive devices [2T1 [22| were fabricated via non-scalable pro¬ 
cesses such as direct laser writing, hampering large scale 
production. 


Recent work has shown that hard polyhedra can self- 
assemble into a great diversity of crystals, liquid crys¬ 
tals, and quasicrystals via entropy maximization [la- 
[35] . a mechanism that can be understood as resulting 
from emergent directional entropic forces (DEEs) be¬ 
tween neighboring particles due to crowding [25l |29l |30| . 
Even chiral crystals have been assembled. The /3-Mn 
crystal structure with chiral space group P4i32 [Fig.[^], 
earlier reported in spheres interacting via an isotropic 
oscillatory pair potential j3T|, self-assembles from achi¬ 
ral dodecahedra shapes m- However, as expected, the 
handedness of the crystal could not be controlled and 
both left- and right-handed crystals were observed with 
equal probability. Given the ability of entropic inter¬ 
actions to imitate the plethora of complex structures 
achievable from atomistic interactions, the absence of ho¬ 
mochiral entropic crystals is both practically discomfort¬ 
ing and theoretically puzzling. 

In this Letter, we investigate whether entropic forces 
alone can be sufficient for the growth of homochiral crys¬ 
tals. By using a recently developed approach for shape 
design [32], we demonstrate that chiral building blocks 
can indeed grow three-dimensional chiral crystals of cho¬ 
sen handedness provided a moderate chiral entropic force 
between neighboring particles can be achieved. We also 
use this method to show how chiral particles can be 
reversibly reconfigured to enforce a homochiral crystal 
when the entropic force strength is tuned by either shape 
‘sharpening’ (or unrounding) or by an increase in deple¬ 
tion concentration. Our results can be rationalized in 
terms of features in particle geometry and suggest design 
strategies currently accessible to laboratory techniques 
that can be used to assemble homochiral crystals. 
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FIG. 1. Assembly of the /3-Mn crystal structure from its 
VPs. (a) The unit cell with 20 particles, (b) The VP shown 
on the left (orange) is elongated with 14 faces and 22 vertices. 
The VP shown on the right (red) resembles a dodecahedron 
with 12 irregular pentagonal faces and 20 vertices, (c) We 
observe the assembly of both left-handed and right-handed 
variants, with a bias for left-handed, the handedness dictated 
by the VPs. (d) A snapshot of the assembled crystal shown 
in projection. The bond orientational order diagram (inset) 
allows identifying the crystal structure. The crystal is degen¬ 
erate, which means either of the two VPs can be found at a 
given lattice point with comparable probability. 


Voronoi particle approach .—The forces controlling the 
self-assembly of non spherical hard particles are emer¬ 
gent and become more directional as particles become 
more aspherical |25[ EH [33]. In fact, by identifying 
shapes capable of maximizing DBFs in a specific crys¬ 
tal, entropically patchy particles can be inversely de¬ 
signed EHl 1301131] ■ One approach is to construct Voronoi 
particles (VPs) [32]. VPs are polyhedra whose shapes 
are identical to the Voronoi cell of the target crystal and 
whose assembly at infinite pressure is, by construction 
and in the absence of degeneracy, guaranteed thermody¬ 
namically to be the target structure. 

We apply the VP approach to the case of /3-Mn in an 
attempt to select the handedness of the crystal. From 
the perfect crystal structure [Fig. Eli we construct its 
two geometrically distinct chiral VPs [Fig. EM , one for 
each set of Wyckoff positions. Although our prior work 
predicts that particles with high sphericity as measured 
by their isoperimetric quotient (IQ) self-assemble into a 
rotator crystal with weak face-to-face contacts between 
neighboring particles [27], /3-Mn VPs have IQ = (0.742, 
0.756), values at the boundaries between crystals and ro¬ 
tator crystals, such that no clear prediction can be made. 
We employ Monte Carlo simulations as in [23] to self- 
assemble the VPs from a disordered fluid state. Systems 
with 2048 particles were equilibrated in the NVT en¬ 


FIG. 2. Assembly of the simple chiral cubic (SGG) crystal 
structure from its VPs. (a) A 2x2 patch of the unit cell with 4 
particles. The bonds indicate chiral helices with pitch length 
four, (b) The VP has point group Us and six symmetry- 
equivalent big faces, (c) The assembled crystal is homochiral 
with handedness determined by the building block, (d) A 
snapshot of the assembled crystal structure shown in projec¬ 
tion. 

semble at packing fractions in the range 0.5 < (j) < 0.7. 
Previous studies of hard particle self-assembly suggest 
this range to be sufhciently high to observe crystallization 
while avoiding kinetic arrest. Structure identification was 
conducted by calculating diffraction patterns and bond 
order diagrams [27] . 

For packing fractions </> > 0.58 we observe that the sys¬ 
tem crystallizes into the target crystal. Just like in sys¬ 
tems of hard dodecahedra [27] and Zetterling spheres [35] , 
we also notice the assembly of y-brass and body-centered 
cubic (BCC) phases competing with the formation of /3- 
Mn. In fact, out of 80 simulations run at these packing 
fractions, only 9 were found to crystallize into /3-Mn, out 
of which 7 assembled into a crystal with the same hand¬ 
edness as the building blocks [Fig. B- Despite this sta¬ 
tistical handshake between building block and bulk crys¬ 
tal chirality, the observation of crystals with unexpected 
chirality suggests that the overall handedness cannot be 
guaranteed a priori and a stronger face-to-face interac¬ 
tion is needed. 

It is worth noticing that we also observe occupational 
degeneracy between the two different VPs [Fig. B]- This 
observation reveals that despite the significant geomet¬ 
ric difference between VPs, their free-energy incentive 
to remain in their respective Wyckoff positions does not 
surpass the gain in mixing entropy |32| [36] . To under¬ 
stand this result, we refer back to the work of Schultz 
et al. [32] where it was demonstrated that the Voronoi 
approach is generally successful if two conditions are sat¬ 
isfied: (i) The unit cell of the target crystal is small, and 
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(ii) the VPs are sufRciently symmetric and of only inter¬ 
mediate sphericity to generate sufficiently strong DBFs. 
If these prerequisites are not satisfied, rotator phases and 
degenerate crystals become possible, as found here. This 
suggests the sphericity of the ,5-Mn VPs might be too 
high to warrant pronounced face-to-face contact, and the 
chiral selectivity is not strong enough for enantioselective 
growth. 

To improve the assembly probability of a chiral crystal 
and to obtain a homochiral crystal, we then hypothe¬ 
size that a more facetted polyhedron, i.e. one leading to 
stronger DBFs, is needed to imprint the chirality of the 
building block onto the assembled structure. As a test, 
we consider a simpler target crystal than /3-Mn. In a 
recent investigation [37] it was discovered that a system 
of spheres interacting via a simple isotropic oscillatory 
pair-potential assembles a cubic chiral crystal with Pear¬ 
son symbol cP4 and the same chiral space group P4i32 as 
/3-Mn. Interestingly, a search through the International 
Tables of Crystallography [38] reveals it to be the only 
crystal up to symmetry equivalence with a chiral cubic 
space group and four or less particles in the unit cell. We 
therefore refer to this crystal as the simple chiral cubic 
(SCC) crystal. Although proposed theoretically in high- 
pressure lithium [39] and semiregular nets [40], SCC has 
not been self-assembled experimentally in natural or ar¬ 
tificial systems. Its low coordination number CN = 6 
together with the low sphericity IQ = 0.564 of its VP 
[Fig. makes SCC a good candidate for the entropy- 
driven assembly of a homochiral crystal. Indeed, we ob¬ 
serve the robust self-assembly of chiral crystals with a 
unique handedness [Fig. in simulations [Fig. [^] at 
packing fractions (f) > 0.60. Between these simulations 
and those of rounded shapes (see below), a total of 92 
simulations assembled into the SCC crystal, all of them 
with handedness dictated by the VP. This result confirms 
that large faceting and low coordination number can lead 
to handedness selection in chiral crystals. A quantitative 
analysis of the strength of DBFs required for such biasing 
will be investigated in the last section of this letter. 

Switching chirality .—Can we systematically modify 
VP shape, and consequently the DBFs, to turn chirality 
off and on at will? Recently it was shown that rounding 
(unrounding) polyhedrally shaped nanoparticles through 
the adsorption (desorption) of ligands can viably switch 
between two different crystals m- To test if this round¬ 
ing approach can be used for chirality control, we modify 
the chiral VPs of the SCC crystal by means of vertex and 
edge rounding [Figs.|^,b]. Rounding is achieved in simu¬ 
lation by using the Minkowski sum of a polyhedron with 
a sphere. Our results show that SCC crystals with the 
same handedness as the building blocks are assembled as 
the particles become increasingly rounded, but only up 
to IQ = 0.743 [22] • Above that value we observe an achi¬ 
ral BCC rotator crystal [Fig. [^]. That we do not find a 
gradual weakening of the preference for homochirality as 



FIG. 3. Reconfiguration between an achiral crystal and a 
homochiral crystal. We compare the assembly of VPs of the 
SCC crystal (a) using their original polyhedral shape, (b) with 
rounded edges and corners, and (c) additionally in the pres¬ 
ence of depletants. (d) Rounding weakens DBFs and leads 
to the assembly of a BCC rotator crystal, (e) The effect of 
rounding can be reversed by depletion, which recovers SCC. 
Depletants were simulated implicitly (see text for details). 
They are shown as a guide in (c) and omitted in (e). 


for /3-Mn, even when the particle sphericity is compara¬ 
ble to that of /3-Mn VPs, can be understood by noticing 
that the low coordination number in SCC precludes the 
formation of a rotator SCC phase. More about this trend 
will be discussed in the next section. 

Another approach to controlling the strength of DBFs 
for colloidal particles is via the addition of depletants [TS] 
12^ . Non-adsorbing macromolecules or smaller colloids 
increase entropic forces between particles that are in close 
contact because it minimizes the excluded volume not 
accessible to the depletants [23]. Recently, depletion has 
been used to tune interparticle distances in colloidal crys¬ 
tals [33] 122] ■ To test whether depletants can reverse the 
trend observed with rounding and re-assemble the ho¬ 
mochiral SCC crystal, we use a recently developed grand 
canonical Monte Carlo algorithm that connects the hard 
colloid system to an external reservoir of penetrable hard 
sphere depletants [21]. The method accounts for deple¬ 
tion effects by sampling how a colloid trial move affects 
the free volume available to the depletants. We combine 
colloids obtained from rounding SCC VPs with small de¬ 
pletants of radius ^ O.lcr [Fig. B] , where a is the circum- 
sphere radius of the colloids. As our simulations demon¬ 
strate [Fig. [^], self-assembly in this system yields the 
original homochiral SCC crystal structure at depletant 
concentrations of > 20% depletants in the reservoir. This 
means the directionality of the depletion force, which is 
strongest at perfect face-to-face contact and thus highly 
directional, is able to undo the effect of rounding and 
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FIG. 4. Quantifying the chiral strength of the particles used in this study, (a-d) PMFT diagrams for the four systems. 
The diagrams were calculated by binning the frequency of finding two particle centers at a given separation vector, creating a 
histogram whose negative logarithmic value corresponds to the energy units reported. The diagrams are two-dimensional slices 
taken through the first neighbor shell parallel to the largest facet on a shape. They were calculated at the lowest densities 
at which each crystal was observed to self-assemble. Color indicates the PMFT, from OkT (dark blue) to —OkT (dark red), 
(e) Coordination Number vs. Isoperimetric Quotient plot derived from Ref. m including the shapes from (a-d). 


restore homochirality. 

Analysis of chirality strength .—We quantify DEFs for 
the four systems presented above by determining poten¬ 
tial of mean force and torque (PMFT) diagrams 
As is visible from the diagrams, /3-Mn VPs have an en- 
tropic interaction of up to 5 kT when faces are in contact 
[Fig.j^a)] . In contrast, the particles assembling homochi- 
ral crystals show stronger interaction of up to 6 kT when 
big faces are in contact (central arrows in the diagrams) 
and up to 4 kT when secondary faces are in contact (satel¬ 
lite arrows) [Fig. ib ,d)]. The rounded particle that as¬ 
sembles only BCC have a weaker interaction of about 
AkT [Fig.|^c)]. In these examples, higher values for the 
PMFT help promote chirality by introducing directional¬ 
ity in the effective entropic coupling between neighbors. 
Geometrical predictions from earlier work 122! are com¬ 
bined with new data in Fig. me). They confirm once 
more that the particles in Fig. ma,c) are too spherical to 
assemble a non-rotator phase, and thus too spherical to 
grow a homochiral crystal. 

In summary, our findings suggest heuristics for the de¬ 
sign of colloidal particles to self-assemble chiral entropic 
crystals with controlled handedness. While a chiral space 
group can occur even for achiral building blocks, mirror- 
symmetry breaking of the colloid is a prerequisite to cre¬ 
ate an enantiomeric imbalance. Sufficiently strong chiral 
forces are necessary for homochirality. Such forces can be 
achieved via particle shape modifications alone, assum¬ 
ing the particle is not too spherical. They can also be 
amplified by depletion effects. Our results demonstrate 
a reconfigurable path for toggling between achiral and 
homochiral crystals via control of an external parameter, 
pushing the envelope towards the design of a new class 
of colloidal materials. 
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